Quantum impurity on the surface of a topological insulator 
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It is shown that the Hamiltonian for a quantum magnetic impurity on the surface of a topological insulator can 
be mapped to the conventional pseudo-gap Anderson impurity model, albeit the combinations of the continuum 
states which hybridize with the impurity have more complex structure in the reciprocal and spin space. If 
the Fermi level is away from the Dirac point, the impurity is predicted to be fully screened at low enough 
temperatures, i.e., there are no residual degrees of freedom. 
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PACS numbers: 73.20.Hb, 72.10.Fk, 72.15.Qm 

Due to the spin-orbit coupling, insulators in the class of 
"topological insulators" have non-trivial topological proper- 
ties of the bulk band structure, which results in the pres- 
ence of peculiar metallic states on their surfaces^. Re- 
cently studied materials such as Bi2Se3 and Bi2Te3 show 
topological-insulator behavior even at high temperatures and 
have the simplest allowed surface states: a single nonde- 
generate Dirac cone^— . These surface states have a char- 
acteristic relation between the momentum and spin, which 
makes the surface electrons insensitive to scattering by im- 
purities since a back-scattering event would require the elec- 
tron spin to be flipped, yet this is prohibited by the time- 
reversal symmetry^—. Recent Fourier-transform scanning 
tunneling spectroscopy (FT-STS) experiments have confirmed 
this predictio n 7 ' 9 ' 10 . Such topological insulator materials may 
find applications in spintronics 5 and quantum computing^. 

The surface-state electrons on topological insulators have 
definite chirality and for the non-degenerate states near the 
Dirac point one has (s(— k)) = — (s(k)) 2 ' ' , Such 
surface states are expected to be perturbed differently by 
non-magnetic and magnetic impurities. The time-reversal- 
symmetry breaking perturbations allow scattering between the 
time-reversed states |k,t) and | — k, The case of clas- 
sical spins (which clearly break the time-reversal symmetry) 
has been explored in Ref. 113[ where it was predicted that the 
impurity opens up a local gap and suppresses the local density 
of states. Scattering on static spin-dependent impurity poten- 
tial was also addressed in Refs. Qj] and EH where character- 
istic interference features were predicted to be observable by 
spin-polarized FT-STS. The description of a magnetic impu- 
rity in terms of classical/static spin is, however, a rather crude 
approximation; magnetic impurities are dynamic objects and 
should be described using quantum impurity models^ 5 -. For 
impurities described, for example, by the Anderson impurity 
modeU£, the local moment of a quantum impurity may be 
compensated at low temperatures by the Kondo effec t 15 ' 17 " — . 
Kondo screened impurities are effectively non-magnetic, thus 
the extent of backscattering on magnetic adsorbates will crit- 
ically depend on the degree of screening of the local moment 
by the surface-state electrons^. Due to the spin-orbit cou- 
pling, there is no SU(2) symmetry in the spin sector alon&» 
therefore it is not immediately clear how effectively the impu- 
rity spin can be screened by the chiral (or helical) surface-state 
electrons. In Ref.[2jl this problem has been approached using 
the Varma-Yafet variational Ansatz. Only partial screening 



(one third) of the local moment was found and 1/T Curie- 
type magnetic susceptibility behavior was predicted at low 
temperatures. In variational methods, there are no a-priori 
guarantees that the correlation functions are correct. In ad- 
dition, the variational Ansatz which projects out the double 
occupancy of the impurity might itself be problematic. For 
this reason, the deduction of the one-third screening from the 
spin-spin correlation functions is uncertain. In this work, this 
problem is reconsidered by a different approach, showing that 
the effective quantum impurity model actually takes a very 
simple form of the conventional single-channel Anderson im- 
purity model, even though the nature of the hybridizing com- 
binations of the conduction-band states is indeed non-trivial. 
For the generic situation with the Dirac point away from the 
Fermi level, the local moment will be screened below the 
Kondo temperature Tk, and the impurity would become ef- 
fectively non-magnetic. No anomalies are expected in the 
low-temperature thermodynamics in this case. Since the ex- 
perimental studies of magnetically doped topological insula- 
tors are already under way22, it is important to arrive at better 
understanding of the Kondo effect in the case of chiral surface 
state electrons. 

In the simplest case, the surface-state electrons on the 
surface of a topological insulator may be described by the 
Hamiltonian^ 



h a p(k) = hvF (kxo-afj + k y o- v a p\ 



(1) 



y } is a two-dimensional (2D) vector in the 



where k = {k x ,k 
reciprocal space of the surface-state band, a and f3 are spin in 
dexes (j" or I), vf is the Fermi velocity, and <j X:V are the Pauli 
matrices.- 8 The a matrix is proportional to real spinal. We 
write k = fc{cos (p, sin </>} where k 
the energy tk = hvpk. We write 



Ikl, and we introduce 



M af3 {4>) ee (cos^a* +an<f>ay) ap = ( J 



-i<f> 



(2) 



and thus /i Q( g(k) = ekM a p{<f>). The hybridization term is 

-ffhyb = ^ka4<Aka + V^lja , (3) 
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where <v and d are the creation and annihilation operator for 
the impurity electron. In this work we take V-^ a = V^; this is 
an approximation, since the surface states labeled by the spin 
index a are not eigenstates of S z . The impurity Hamiltonian 
is simply 



imp 



(4) 



where n a = d' a d a is the occupancy of the spin a level, e is 
the impurity energy level, and U the on-site electron-electron 
repulsion. 

The reduction of the problem to an effective one- 
dimensional quantum impurity problem is analogous to the 
derivations in Ref. |24l but performed for the 2D case. We go 
from the discrete k to the continuous vectors k in the stan- 
dard way, / 2 „p / d 2 k, where N is the number of 
the surface states, and the components of the wave vector on 
the right-hand side run from — ir to it. The continuum opera- 
tors that correspond to are Oka and they are normalized 
as {aka)«k'a'} = ^(k — k')<y o a /, thus the mapping of the 
operators is -0ka ~> (27r/-\/]V)ak Q . This leads to 



H = y^l d 2 kaj CQ /i a( a(k)a k /3, 
H hyh = j d 2 \t V k dia ka + H.c. 



In the next step we expand the operators a kQ in the az- 
imuthal components. In fact, we allow for a slightly more 
general Ansatz: 



V« V27T 

m,p 



(6) 



where U is some unitary matrix which may depend on the 
azimuthal angle <f>, while m € Z. The operators Cfe TOy g are 
normalized as 



{c/cm/3,4, m ,^,} = <5(fc- fc')<5m,m"5/3^'. 



(7) 



The inverse transformation is 



Choosing U = W in Eq. ©, we obtain a simple diagonal 
form for the band Hamiltonian: 



H o=J2j dke k<J c\ miJ c km<J , 



(10) 



with efrf = efc and = —e k . The density of states is 



P<r(e) 



(2^ 



d 2 k5(e - e kc 



N 1 



(ii) 

where ^t(x) = 6{x) and #^(2;) = 6(— x), with 0(a;) the Heav- 
iside step function. Thus a —\ corresponds to the upper Dirac 
cone, and a =1 to the lower Dirac cone. 

For simplicity, we will at first assume the hopping constants 
V k to be isotropic in the 2D space, i.e., Vk = V k . Such k- 
dependence corresponds, for example, to the d 3z 2 „ r 2 impurity 
orbital. The hybridization term then transforms as 



ffhyb = Yl [ ^dkV k 



dlckmfj + H.c. 

(12) 



(5) Noting that for integer m 



1 1 
V52^ Jo 
we obtain 
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d<f> e 



i(m±l/2)4> 



1 1 



4/ 



y/2 2?r 2m ± 1 



= lt (13) 



Hhyb — 



N 



+ 



J VkdkV k 

d^Ckrnf} + H.C. 



(14) 



a/3 



Finally, we introduce the energy representation by defining 
Cerna- = (de k cr / 1 dk)~ l ' 12 [ c kma , with e = e k(J . The normal- 
ization then becomes {c £mCT , c^, m(T ,} = 6(e - e') 5 m,m'^,<r'- 
Note that e > for a =t and e < for a =\,. We introduce 
an upper energy cutoff for a =f at +£> (the lower limit is 0) 
and a lower energy cutoff for a =1 at —D (the upper limit is 
0). We then obtain 

D r 



The matrix M in Eq. (O may be diagonalized using the 
unitary matrix 



(15) 
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e'^ 

,i4>/2 Q 



-ia y 7r/4 
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X f e -i<t>/2 e -i<t>/2 
a i<)>/2 _ p i<}>/2 



(9) 



For the hybridization term we have 

[ r- 

Vk dk VkCkm/3 

0v 



/27T 



which gives W^MW = ^ 1 = cr 2 . The transforma- 
tion matrix has period 47r, which reflects that a quantized 
Berry's phase of ir is acquired by an electron circling the 
Fermi arc, as characteristic for a "topological metal" surface 
phase where the Fermi arc encloses a single Dirac point^. 
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k f de 



dk 
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(16) 



'2tt J " \d€ kfj 
= j de[pp(e)]^ 2 V(e)c imP , 



V(e)c em p 
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since the density of states is p a (e) = (N/2ir)k e \dk/de\. We 
conclude that 



m,ct 3 $ 



Trn 

-7^ 



c4c em/3 +H.c. 



a/3 



(17) 

The hybridization function is defined as r cr (e) = 
7rpcr(e)|F(e)| 2 . Since the coupling to a continuum of states 
in a quantum impurity model is fully defined by its hybridiza- 
tion function r(e), we replace */ pp(e)V(e) by y/Tp(e)/Tr in 
the following. 

We now introduce the combinations of states 



h f _„ = 



It. 



1m c e 



7m c e 



(18) 



where the normalization factor r is defined as 

-I 1/2 



E 



7n 



= 1 



(19) 



The two sets g and h are (canonical) fermionic operators, 
{g ea , gti-,} = S(e — e')5 aa , and likewise for h, and they are 
orthogonal to each other: {g e <j, hi, ,} = 0. 
We obtain 



H = dee g*g et + dee g\,g ei 

JO J-D 

+ dee fiLhtf + dee hlh £ i, 

JO J-D 

H hyh = J de s/Y(e)/-K r 

x (0^{e)d\h et + 0i(e)d\h cl 
+ e t (e)d\g et - 6 x {e)d\g ei 



(20) 



H.c. 



We introduced T(e) = T^-(e) + r^(e), which is non-zero for 
all e except at e = 0, where there is a linear pseudo-gap in the 
density of states. In Hq we omitted writing those combina- 
tions of the conduction-band states which do not couple to the 
impurity, as they play no role in the following. 

We note that the impurity spin-up orbital only couples to 
the h states, and the impurity spin-down orbital only to the 
g states. Since the g and h states are only defined for the 
combinations e > 0, a =t and e < 0, a =1, we may actually 
drop the spin label and write 



9t 



9 t (e)g et - 6i(e)g eh 

e t {e)h ef + e^h^. 



(21) 



Thus we are dealing with two channels of spin-less electrons 
(defined for all e) with a pseudo-gap density of states. The 



Hamiltonian is then 

,-D 



Ho = I dee (g\g e + h\h e ) 

J-D 



H hyh = de VT(e)/n r [d\h t + d\g^j + H.c. 



(22) 



As a final step, the particle type label h and g can be replaced 
by a (pseudo)spin label, i.e., = h t and b e i = g e . The oper- 
ators b ea are canonical fermionic operators. The final Hamil- 
tonian takes the form of the regular single-impurity Anderson 
model (compare with Ref.l24f): 



i?0 = E / de6 b tcrbea, 
- J-D 



#hyb = E / de VTWA T d\b ta + H.C. 
„ J-D 



(23) 



The end-result of the derivation is thus trivial in its form, 
as the low-energy effective model corresponds to the conven- 
tional Anderson impurity model (with a pseudo-gap in the hy- 
bridization function). Nevertheless, it should be noted that 
the operators b ea correspond to a rather complicated combi- 
nation of the original conduction-band electron states. Re- 
cently, formally the same effective model has been studied in 
the context of Anderson impurities adsorbed on a graphene 
layer—. It is interesting to note that despite the differences 
between the physical content of the Hamiltonian in these 
two different cases (as already observed in Refs. [TH andl2lh 
the effective impurity model in a suitably transformed ba- 
sis is the same. It should also be noted that this reduc- 
tion is only possible starting from an Anderson-type impu- 
rity model. Starting from a Kondo-type impurity model, 

H K = JS ■ {J2kk' a p ^la^a^k'p), we would end up with 
an effective model with an infinite number of channels (index 
to), because there is no equivalent of Eq. ( fT8l to perform a 
separation of the linear combinations of states. If a spin-only 
model is desired, the Schrieffer- Wolff transformation has to 
be performed as the last step of the derivation. 

For orbitals which couple with the continuum via Vk which 
depends on the azimuthal angle, a similar result would en- 
sue. For example, for d xz and d yz we have Vk = Vk cos <j> 
and Vk = Vk sin <fi, respectively, while for d x 2_ y 2 and d xy we 
have Vk = V/ s cos2</> and Vk = Vfcsin20, respectively. It 
easy to see, however, that the angular additional factors cor- 
responds to a mere shift of to by 1 or 2 in Eq. ( fT3b . thus we 
end up with the same effective model, although with different 
combinations of continuum states coupling to the impurity. Of 
course, the Fermi surface/line is not axially symmetric either, 
but rather reflects the symmetry of the lattice (for example, 
due to warping effects resulting from the higher-order spin- 
orbit coupling terms, the Fermi surface in the case of Bi2Te3 
has a hexagram/hexagon shape and becomes circular only in 
the vicinity of the Dirac point 5,8 ). Again, this only affects the 
combinations of states which couple with the impurity orbital, 
however it does not change the results in a qualitative way. 

The properties of the Anderson impurity in pseudo-gap 
Fermi baths are well know n 26 i 27 : for a linear pseudo-gap, 
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p oc \us\, the system flows to a local-moment (LM) fixed point 
where the impurity is effectively decoupled from the conduc- 
tion states and remains unscreened at low temperatures, unless 
the hybridization T is strong enough. The LM fixed point for 
€p = is very resilient to various perturbation; in particu- 
lar, it persists for a large range of parameters away from the 
particle-hole symmetric point 2 - 7 -. With e p away from the Dirac 
point, however, the density of states is finite and the impurity 
is screened in the conventional Kondo effect, thus the system 
flows to the strong-coupling fixed point. In topological insu- 
lators, unlike in graphene, there is no a-priori reason for the 
Fermi level to sit at the Dirac point, thus the impurity will 
typically undergo Kondo screening at some low-enough tem- 
perature Tjc- The impurity contribution to the total entropy 
Simp is much below T g 15,19 ' 24 , which implies that despite 
a complex structure in the spin/reciprocal space, no anoma- 
lies are expected in the low-temperature thermodynamics (in- 
cluding the magnetic susceptibility): the impurity spin is fully 
screened (not one-third screened^). This does not preclude 
a complex structure of the Kondo cloud, which might indeed 
exhibit non-trivial spatial and spin dependence. It is, how- 
ever, unclear whether these features could ever be directly ob- 
served, since the Kondo cloud is rather elusive. 

An Anderson-model-type magnetic impurity on the sur- 
face of a topological insulator may be fully Kondo screened. 



For temperature much below the Kondo temperature Tk, the 
backscattering of the surface state electrons will be prohib- 
ited, as in the case of non-magnetic adsorbates. No open- 
ing of the gap in the surface-state band is expected in this 
regime. At temperatures above Tk, however, the spin-flip 
scattering events will connect the points on the opposite sides 
of the Fermi surfaces. The possible opening of the gap 
might be, however, masked be thermal smearing effects. On 
the other hand, for ferromagnetically Kondo coupled impu- 
rities or for large-spin impurities with easy-axis magnetic 
anisotropy, Kondo screening is non-effective and strong scat- 
tering effects may indeed be expected. Further work should 
explore the behavior of multiorbital impurity models and take 
into account the full spin and orbital structure of the hybridiza- 
tion parameters Vka. m p, where m is the impurity orbital in- 
dex, while a and j3 are the spin indexes of surface and impu- 
rity states, respectively. 
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